Transition to a Layered Phase-Decoupled Superconducting State in La 2 _ :r Ba x Cu04 
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We present measurements of the resistivity parallel and perpendicular to the Cu02 layers in single 
crystals of I^-^Ba^CuCU (LBCO) with x — 0.095 as a function of perpendicular magnetic fields 
up to 35 T. At temperatures below the zero-field superconducting transition of 31.5 K, there is a 
regime of field strength for which superconducting phase coherence between the layers is destroyed, 
yet effective coherence is retained within the layers. Evidence that this is an actual state of matter is 
provided by a well-defined loss of order, indicated by the field dependence of the in-plane resistivity. 
The behavior of this phase is qualitatively identical to that observed previously for LBCO with 
x = 1/8 in zero field, in which frustration of the interlayer Josephson coupling is associated with 
spin and charge stripe order. 



Unusual superconducting states have been observed [JJ 
[2] in association with the 1/8 anomaly in La2_ K Ba K Cu04 
(LBCO) OHj. In this earlier work, the three-dimensional 
(3D) superconducting transition temperature, T c , is de- 
pressed below 5 K; nevertheless, cooling below ~ 40 K 
yields an order of magnitude drop in the resistivity par- 
allel to the Cu02 planes, pii . On further cooling, pu de- 
creases in a fashion consistent with that expected for a 
two-dimensional (2D) superconducting vortex liquid [5], 
while the resistivity perpendicular to the layers, p±, re- 
mains large; this apparent layered-vortex-liquid (LVL) 
state is accompanied by anisotropic diamagnetism. Fi- 
nally, pii extrapolates to zero at ~ 16 K, while p± remains 
finite. Regardless of whether the individual layers de- 
velop true superconducting order, one would expect that 
inevitable Josephson coupling should have driven 3D su- 
perconducting order near 40 K where the first strong 2D 
superconducting correlations appeared [6]. The strong 
suppression of 3D order indicates that the interlayer 
Josephson coupling is frustrated for some reason [THS]- 

While the frustrated Josephson coupling is surprising, 
the occurrence of superconducting order in a system of 
decoupled layers is compatible, in principle, with existing 
theory [10], at least in the absence of an applied mag- 
netic field. What is more surprising is that both the LVL 
and the layered phase-decoupled (LPD) superconducting 
(SC) state survived in a perpendicular magnetic field, 
/xo-ffj_, of at least 3 T. Such a field should easily destroy 
the order in individual layers |llj . Recently, transport 
studies of LBCO with x = 0.095, a good bulk supercon- 
ductor, indicated a field-induced LPD-SC state [121 U5] : 
however, with the limited range of magnetic field used 
it was not possible to distinguish the thermal variation 
of the resistivity from the exponential behavior expected 
for thermally-excited vortex motion [Til HI] • 

In this letter, we report transport studies of LBCO 



x = 0.095 in poH± up to 35 T. The data provide evi- 
dence for a high-field thermal crossover from the normal 
state to an LVL phase, followed at lower temperature by 
a transition to the LPD-SC. Our results are summarized 
in Fig. [TJ In zero field, the sample has T c q = 31.5 K in 
measurements of both p± and py , but the observed transi- 
tions separate at finite field. The onset of finite p± occurs 
at H^r , while, as we will explain, phase-decoupling of the 
layers happens at the slightly higher field Hq. In con- 
trast, we observe evidence for a phase transition to finite 
p\\ at He, with H'j S> H^r. The temperature dependence 
of p\\ above He is similar to that predicted [51 [TS] for 
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FIG. 1. (color online) Phase diagrams in H±-T space obtained 
from measurements of (a) p± and (b) p«. In (a), triangles in- 
dicate the onset of finite p± at H^r ; squares denote Hq corre- 
sponding to the interlayer phase-decoupling crossover. In (b), 
circles indicate onset of finite p\\ at Hi ; vertical solid line cor- 
responds to T 2D , the crossover from the layered vortex liquid 
phase to the normal state. In each panel, the shaded contours 
correspond to the resistivity normalized to an extrapolation 
of the normal-state behavior obtained at the maximum field. 
(Online: p/p n < 10 -3 is cyan.) 
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FIG. 2. (color online) Measurements of (a) pj_ and (b) pu as 
a function of H±, obtained at various fixed temperatures as 
listed in each panel. 



a two-dimensional (2D) superconductor above a vortex- 
antivortex unbinding transition [IH[T7], Ibkt- At T^ D 
there is a crossover to the normal state. The state be- 
tween Hq and He we identify as LPD-SC, while the state 
between He and T^ D is identified as the LVL. 

The crystals, grown by the traveling-solvent floating- 
zone method, have been characterized in several previous 
studies [U Q21 Q21 [18] . The present experiments were 
performed in the 35-T dc magnet at the National High 
Magnetic Field Laboratory (NHMFL). The crystals for 
the measurements of p± and py are the same as those 
used in a previous transport study [12] . and the contact 
configurations are described there. The resistance was 
measured using an ac resistance bridge with an excitation 
current of 1 mA. All measurements were done after field- 
cooling from above T c q. The resulting scans of resistivity 
vs. H±_ obtained for a range of temperatures are shown 
in Fig. [J 

Let us first consider p±(H±). For T < 30 K, increas- 
ing H± initially causes p± to become finite, followed by 
rapid growth and eventual saturation, followed by a grad- 



ual decrease beyond the maximum. Such behavior has 
been studied previously, especially in F^S^CaC^Os+a 
pi 120] and Bi 2 Sr2_ x La a; Cu06+ ( 5 [21]. The rise of p± 
with increasing H± is due to suppression of the conduc- 
tion channel associated with interlayer pair tunneling; 
on crossing the maximum, single-particle transport dom- 
inates [Jj5|. The region of negative magnetoresistance 
at high field has been attributed to the impact of H± 
on the pseudogap [20] ; reducing the antinodal gap in- 
creases the density of normal carriers that can move be- 
tween planes. Parallels have also been drawn with the 
field-tuned superconductor-insulator transition observed 
in disordered thin-films of various metals [22 24 . By 
this latter analogy, the resistive transition in p± can be 
viewed as a transition to a Cooper pair insulator phase 
at high /io-ffj_. In our case, the Cooper pairs are localized 
along the c axis, becoming restricted to the Cu02 layers. 

Is it possible to identify a threshold for insulating be- 
havior along the c axis? It has been argued by several 
groups that the rise in p± can be understood in terms of 
phase fluctuations in the interlayer Josephson junctions 
due to thermal noise [2"BTt28| . In this interpretation, the 
relevant quantity is the extensive resistance per Joseph- 
son junction. Hettinger et al. [3B] demonstrated em- 
pirically that the effective junction area corresponds to 
A = $o/(-B_l +B ), where $o is the flux quantum and B 
is a parameter. We have shown in [T^] that this approach 
gives a good description of the evolution of p±(T, H±) in 
our sample for T < T c q with B± m poH± and Bq = 2.2 T. 
The effective junction resistance is then R± = p±s/A, 
where s is the interlayer spacing (6.6 A). 

According to Halperin et al. [29], the criterion for the 
Josephson junctions to become effectively insulating is 
R± > Rq, where Rq = /i/(4e 2 ) = 6.45 kfl is the quan- 
tum of resistance for Cooper pairs. We define Hq as the 
field at which Rj_ = Rq . In Fig. [3j we show that Rj_ /Rq 
scales as [(H±_ - H^)/(H Q - Hj)] a(T \ with the T de- 
pendence of the exponent a displayed in the inset; the T 
dependences of H^r and Hq are shown in Fig.JlJa). The 
scaling is motivated by a calculation from Konik [30| for 
p± in a model of weak Josephson coupling between 2D 
layers; he predicts a K = |[1 + (T c0 - T)/T c0 ]. To par- 
tially compensate for changes in the field dependence on 
converting from p± to R±, the dotted line in the inset 
corresponds to cvk + 0.5. 

In Fig. g] we compare the temperature dependence of 
p± and pti in high fields. For p±, the field appears to 
shift the superconducting transition to low temperature. 
In contrast, pu shows a substantial drop near 30 K even 
in the highest field, and it continues towards zero on fur- 
ther cooling. To emphasize this behavior, the inset of 
Fig.[4|b) shows the ratio of finite differences Apy/AT, as 
a function of field and temperature. For poH± ^ IT, one 
can see that the maximum of this approximate derivative 
occurs at 30 ± 2 K, which we identify as T C 2D , the onset 
of mean-field 2D-like superconductivity [3TJ. The finite 
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FIG. 3. (color online) Plot of R±/R Q vs. [(H± - H^r)/(H Q - 
H^)] a ^ T \ as described in the text, for data at temperatures 
from 1.6 to 27 K. The dashed straight line is for reference. The 
T dependence of the exponent a is plotted in the inset. Dotted 
line represents «k + 0.5, as described in the text; solid line 
corresponds to a = l.l + 0.75t + 4.5t 2 , with t = (T c0 -T)/T c0 . 



resistivity at lower temperature indicates a lack of phase 
order, as we discuss next. 

In zero field, a 2D superconductor has zero resistiv- 
ity only below the Bcrczinskii-Kostcrlitz-Thouless tran- 
sition, Tbkt, as above that temperature depairing of 
vortex-antivortex excitations destroys the phase coher- 
ence [IHl HZ]- The resistivity, which is proportional to 
the product of the density of free vortices and the vortex 
mobility, is predicted [SJ [T5] to have the form 



P\\/P\\ 



-b/T 



(1) 



where p» n is the normal-state resistivity, r = T/Tbkt — 
1, a and b are constants of order one, and 7 = 0.5. In 
the mixed state, the applied H± will increase the den- 
sity of free vortices, and hence will increase the resistiv- 
ity [3H [33J . In a study of 2D-like superconductivity in 
LBCO x = 1/8 PQ, it was found empirically that P\\{T) 
in a field could still be described by Eq. 0, provided 
that one takes account of the reduction of Tbkt by the 
field. Theoretically, one does not expect Tbkt to remain 
finite in the mixed state of a 2D superconductor [IT] : 
however, our system is never truly 2D. We have a stack 
of 2D layers, and the electromagnetic interaction between 
pancake vortices in neighboring layers provides some in- 
terlayer coupling, allowing a BKT-like transition in the 
zero-field limit 10 . 

To apply Eq. Q, we need an estimate for py )n . We 
choose to use a linear extrapolation of p^ fit to data at 
T > 50 K for pqH± = 35 T. Using this, we obtain the 
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FIG. 4. (color online) Measurements of (a) p± and (b) pu as 
a function of temperature for poH± = 0, 20, and 35 T. The 
dashed lines in (a) are guides to the eye. The solid lines in 
(b) are calculations of 2D flux-flow resistivity, as discussed in 
the text. Inset of (b) shows Apy/AT as a function of H± and 
T. 



curves through the 20 and 35-T data in Fig. gb), with 
parameters Tbkt = 16.5 and 12.5 K, respectively, a = 2.5 
and b = 1.5; here and in the following we set 7 = 0.6 to 
improve the agreement with the data. We can use the 
same formula to compare with p\\(H±)/p\\ n of Fig 
provided that we insert the field dependence throug 
T B kt(#l); the latter corresponds to ffj (T) of Fig.^b). 
We get an initial estimate Hc l (T) from the field at which 
p\\/p\\ t n falls below 10 -3 , yielding the (black) squares in 
Fig. [5] To account for the finite cutoff, we make the 
adjustment p H 1 } = p, Hf - d/T with d = 44.8 TK. 
Inverting, we fit Tbkt with a cubic poynomial in H±. 
The curves calculated with Eq. (JTJ) are indicated by the 
solid (black) lines in Fig. [5j where we have used the same 
fixed values of a, b, and 7. Although there is no theo- 
retical justification for it, we find that Eq. (JlJ provides a 
consistent fit to the data with few parameters. 

To review, we see in Fig.|4]that cooling from the normal 
state in a large H± results in a drop in pn at T^ D w 30 K, 
while the temperature dependence of p± has insulating 
character. The drop in p» must be due to superconduct- 
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FIG. 5. (color online) Plot of p^(H±) normalized to approx- 
imate normal-state values (see text). Squares indicate fields 
at which data cross I0" 3 . Solid lines are calculations using 
Eq. |T]), as discussed in the text. 

ing correlations, which are only apparent for transport 
parallel to the CuC>2 layers. Thus, it seems reasonable to 
interpret Tj? D as the onset of mean-field superconductiv- 
ity within Josephson-decoupled 2D layers, just as occurs 
at 40 K (in zero field) for LBCO 1=1/8 0. With further 
cooling, or on reducing the field, the resistivity decreases 
in a fashion that is qualitatively consistent with expecta- 
tions for vortex-liquid behavior in the layers. The sharp 
drop in on reducing H±, shown in Fig.[5j strongly sug- 
gests a transition from a vortex liquid to a vortex solid 
state [34] . 

A vortex-liquid state has previously been proposed to 
explain features of the pseudogap phase at T > T c q 
[551 I5B"] . Along with other recent work J37J 13H] , our re- 
sults suggest that such a scenario is overly optimistic. 
The rise of py on warming through T C 2D indicates a loss 
of uniform superconducting correlations in the normal 
state. T c q is slightly larger than Tj? D , suggesting that, 
upon cooling in zero field, 3D order develops just before 
the individual layers would become superconducting in 
the absence of Josephson coupling, as others have ob- 
served in Bi 2 Sr 2 CaCu20 8+1 5 HMD- 
Observation of the LVL (and the LPD-SC) is only 
possible when the interlayer Josephson coupling is frus- 
trated. We note that one proposed origin of the frus- 
tration is the development of pair-density-wave (PDW) 
order in association with charge and spin stripe order 
|H1 H2] • Neutron and x-ray diffraction measurements 
on LBCO x = 0.095 have demonstrated that both charge 
and spin stripe order are enhanced by H± for T < T 2D 
[12]; however, it should be noted that, while the occur- 
rence of PDW order would explain the loss of 3D su- 
perconducting order, it would not, by itself, explain the 
apparent stability of the LPD-SC state in large H±. 



The modulated pair wave function of the PDW state 
provides a way for the superconductivity to coexist with 
local antiferromagnetic order: the superconducting wave 
function has zeros at the extrema of the spin density 
wave, and vice versa. This is consistent with experi- 
mental evidence that long-range commensurate antifer- 
romagnetic order and superconductivity do not coexist 
in LBCO or La2- a; Sr :r Cu04 [43U45] , At the same time, 
the modulation makes the PDW state quite sensitive to 
disorder, consistent with the strongly depressed bulk T c 
in LBCO x — 1/8. For x = 0.095, the superconduct- 
ing order develops at much higher T c q and T^ D is vir- 
tually independent of H± even for (j,qH± as high as 35 
T. The stability indicates that there are at least quan- 
titative differences from x = 1/8. Nevertheless, it will 
be shown elsewhere [IB] that, even in zero field, the low- 
energy incommensurate spin fluctuations are gapless and 
that there is no sign of a conventional spin resonance. Is 
there any solution other than a spatially-modulated su- 
perconducting state that could coexist in this material 
with gapless, spatially-modulated spin fluctuations? 
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